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ABSTRACT 



The temperature dependence of diffusion length for 
PU-BE neutrons in distilled water poisoned with boric 
acid was studied using concentrations of 1.0, 8.0, 16.0, 
and 30.0 grams per liter. Runs at each concentration were 
made at temperatures of 95, 131, and 167° F. The trend of 
temperature dependence was clearly established. A theo- 
retical temperature dependence of k^, where k = 1/L is- 
proposed using a Maxwellian distribution and found to be 
in good agreement with the data reported in Reference [2]. 
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I. 



INTRODUCTION 



The dependence o£ thermal neutron diffusion length 
upon poison concentration in water has been the subject of 
several experiments, the most extensive of which was done 
by Goddard and Johnson [2]. Using boron, a strong 1/v 
absorber- - its cross section varies inversely with neutron 
velocity- - these investigators carried out experiments over 
--a wide range of boron concentrations. 

The object of this study was to determine diffusion 
length as a function of boron concentration and tempera- 
ture. An attempt was made to represent both the absorp- 
tion and the temperature dependence of diffusion length in 
a polynomial approximation. This experiment was conducted 
using a plutonium-beryllium source. 
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II. THEORY 



A. ONE GROUP THERMAL NEUTRON DIEFUSION 

For a point source of thermal neutrons in an infinite 
hydrogeneous medium, the one grou]) thermal diffusion equa- 
tion may be written in spherical coordinates as 



Using 



r^ dr 



the source 



lim r^ 
r->“ 



d 

. 2 i. 

dr 



condition 



J (r) 



S_ 

4tt 




( 1 ) 



where , 



J 



-D 



dr 



the solution to the diffusion equation is, 

*^T 4irDr (2) 

Where S is the source strength in neutrons/sec. 

The determination of the diffusion length v^ould now 
be possible using the above solution if the neutron source 
were, in fact, thermal. Hoivever, the PU-BE neutron source 
emits neutrons with an average kinetic energy of approxi- 
mately 4.2 MEV. The actual neutron spectrum extends from 
very low energies to about 10.6 MEV. In order to use 
the one group thermal model, it is necessary to account 
for the PU-BE neutron energy spectrum. This is accomplished 
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t 



by considering the neutrons as being born in two groups ; 
a fast group with an energy range above .025 ev up to 10.6 
MEV and a thermal group with energy .025 ev. The next 
simplifying assumption is to assume all the neutrons born 
to the fast group are born at the average group energy. 

The actual energy distribution is shown in Figure 1, hence 
there would be many different migration areas [6]. It is 
assumed that the fast neutrons are born at an average 
energy and therefore, there will be an average migration 
area. The average energy of the PU-BE neutrons is 4.2 MEV. 

Now having lumped the fast neutrons into a single 
energy group, the source term in tlie diffision equation can 
be derived. This term must represent those source neutrons 
that are born fast and have slowed down to thermal energy. 
The one group diffusion equation written with the source 
term has the following form; 



B. SOURCE TERM DEVELOPMENT 

For an isotropic point source in an infinite medium. 



hydrogeneous media a fast neutron may loose all energy 
after one collision, the density of first scattering 
collisions is given by [5] . 




(2A) 



the source releases S neutrons per unit time. Since in 

o ^ 




Z^(EQ)e'^s^^0^ 



(3) 
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For an infinite hydrogen medium, the slowing down 
density at any energy is constant and equal to the source 
density [4] 



S(r) = q (r,E) 



(4) 



Therefore, for a point source emitting fast neutrons at 



E^, the number of neutrons entering the thermal region is 



given by 



S(r) = q (r, E^^j) 



S S (E ) e 

O S ^ O-^ 

4irr^ 



-£sCEo)r 



substituting into equation (2A) results in 



1_ ^ + ^o 

r^ cIF dr L.2 4TrDr^ 



(E^)e 



-IsCV^ 



= 0 

( 6 ) 



or 



1 _ _d_ 
r ^ dr 



d4>™ 
.2 L 



dr 



1 + C e'^/^ = 0 



(7) 



where , 



C = 



S E(E ) 
0^0^ 

47rD 



and 



l/b = 



C. SOLUTION OF DIFFUSION EQUATION 

In the development a first collision density is as- 
sumed and is assumed constant at the average energy of 
the fast neutron. The homogeneous solution to this equa- 
tion is: 






H 



A 




r 



+ 



. + r/L 

B 

r 



( 8 ) 
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FIGURE 1. Neutron Spectrum for a PU-BE Source. 
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Since the flux must be finite everywhere , B must be zero. 
Assuming a particular solution of 



P 




(9) 



and substituting in Equation (7) yields 




-r/L 



Dr"=(L^-b2) 



T 



r 



( 10 ) 



This is the asymptotic solution (i.e., the solution for 



r->«>) for Equation (6) . 

If is larger than b, the first term of Equation (10) 
becomes much more significant than the second term at some 
sufficiently large r and the flux can be assumed to decay 
as a simple exponential. For PU-BE sources b is larger 
than , hence the complete expression must be retained 
to correctly describe the neutron diffusion. The several 
unknows associated with this equation; namely, the con- 
stant A, the parameters D,b considerably increase the 
difficulty of determining an accurate diffusion length. 

Even though it has been possible to develop a theoretical 
expression for flux based on a one group model with a 
source term to account for the fast neutrons, this ex- 
pression is not simple and does not lend itself to a 
straightforward determination of diffusion length. How- 
ever, it is still possible to proceed by turning to the 
relationship between diffusion length, age, and migration 
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area. Although this method is far from simple it does 
offer a solution. The next step will be to develop the 
expression for age and migration area. 

D. NEUTRON AGE CALCULATION 

The definition of age can be derived from the reactor- 
criticality equation written in the moment form [6] 

f C-D" 

1+L^B“ ° (2n+l) ! (12) 

For an infinite medium it is possible to neglect all terms 
greater than n = 1. The summation of equation (12) can be 
written as : 

T B^r^ ^ B'^r'* , B®r® , 

I - . + 4- + • • • 

3! 5! 7! 

Now setting the above approximation equal to the critically 
equation in the form 



Koo e 



1+ L^B 



•B=^t 



2 152 



= 1 



(13) 



we get 



-B^t , B^r^ 
e = 1 - 



3! 

when all terms in the series with n>l has been neglected. 
Using a series expansion of the exponential, the expression 
becomes 



1 - B^-t = 1 - 



2T2 



B^r 

3! 



(14) 
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which yields : 



(Eo,E) = 






( 15 ) 



The age relates to the mean-square distance traveled by a 
neutron from its point of birth as a fast meutron to the 
point at which it becomes thermal. 

The mean- square (net vector) distance, r ^ , which a 
neutron travels from its source to where it is absorbed, 
is given by [3] . 

r 

I r^ 4in:^E <i> dr 

Yz = ^ 

fOO 

I 4iTr^ E $ dr 

Jo a (16) 



Substituting 



= 3_il^ 

^ s 



(17) 



from the age-diffusion theory, gives 



-2 2 q (r ,E) dV 

^ / q (r,E) dV 



(18) 



The saturated activity can be written as 



A - V ^ 
s a 



(19) 



where both the macroscopic cross-section and the flux are 
functions of energy. Integrating over the energy range of 
interest and substituting Equation (17) into Equation (19) 
yields : 
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q (E) E^(E) dE 



C E E 

^ s 



( 20 ) 



For the single resonance peak o£ indium it is pos- 
sible to assume: (1) the integral is small except at the 

resonance peak, (2) the slowing-down density is independent 
of energy and (3) the moderating power (C^g) is constant 
over the peak. Equation (20) may now be written as 



This shows that the slowing down density is directly pro- 
portional to the saturated activity. Substituting Equa- 
tion (2) into Equation (18) yields: 




( 21 ) 



r 



2 



r^ Ag (r) dr 



s 



A (r) dr 



( 22 ) 



which for cadium- covered foils becomes 




r" A^p (r) dr 



r^ (r) dr 



(23) 



Finally 




0 



( 24 ) 
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An experimental evaluation of Equation (24) is in- 
cluded in Appendix B. 

E. DETERMINATION OF MIGRATION AREA 

The total mean-square distance a neutron travels from 
its point of birth to its point of interaction, con- 
sists of two parts: (1) the mean-square distance covered 

in slowing down process, r^^ and (2) the mean-square dis- 
tance covered in the process of thermal diffusion, 

The first parameter correlates with the age, and the second 
parameter correlates with thermal diffusion length, [6]. 
The migration area is defined as 



M^ = T + L- 



(25) 



The correlation betiveen r^^ and M^ is given by 



r** 

r^ A (r) dr 

(26) 

The integrals can be obtained in a manner analogous to the 
derivation of Equation (24) for 

F. TEMPERATURE DEPENDENCE APPROXIM/VTIONS 




In order to formulate the temperature dependence of 
diffusion length a series expansion in absorption cross 
section Z of the form suggested by reference [2] was 

3. 0 

utilized . 






ao+ ai 



azY.^ 



ao 



+ 






3 

ao 



( 27 ) 
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where K is the inverse o£ the diffusion length, Z , is 
the sum of the 2200m/sec macroscopic absorption cross sec- 
tions for boron and hydrogen and the a's are coefficients 
to be determined. We shall attempt to include the tempera- 
ture dependence in the expansion defined by Equation (29) . 

A first approximation of temperature dependence was 
derived by assuming a Maxwellian distribution [4] for 
neutron flux. 



where n is the density of neutrons, K is Boltzmann's con- 
stant, m is the meutron mass, and T is the temperature of 
the medium in degrees Kelvin. By setting the derivative 
of Equation (28) equal to zero it can be shown that the 
most probably velocity is given by [4] 



where is the most probable velocity for a Maxwellian 
flux at 293. 6°K. 

Using this relationship it is possible to relate the 
average thermal macroscopic cross section as a function of 
temperature [4]. 




(28) 



V.p = 1.284 Th X 10“ ‘^’^^sec (29) 



from this 




(30) 



16 



f 44^ i i A ■■■ ' " 






i 








V 



( 31 ) 



/ TT 



E = 
a 






o 



S (E ) 
a ^ 



This relationship v/as applied to the polynomial approxi- 
mation with the temperature dependence being applied to the 
coefficients. In order to test the validity of the approxi- 
mation the data reported by reference [2] at 22.3°C was used 
for the initial data at T^. This approximation was found to 
correlate very well \\'ith the data at 65°C from the same 
reference. The relatively large error associated with the 
present' research precludes any conclusion more specific than 
that trend of the data is in the direction of the approxi- 
mation. Figure (7) is a plot of versus E for various 
temperatures using the temperature dependent approximating 
polynomial . 
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III. DESCRIPTION OF THE EXPERIMENT 



Indium foil activation measurements were made in 
heated water poisoned with boric acid (H^BO^) . Fifteen 
indium foils were placed in a plexiglass holder at inter- 
vals of 1/2 inch. The holder was located in the center of 
a 22 X 30 X 23 inch tank made of 1 inch plexiglass. The 
tank contained approximately 45 gals, of distilled water. 

Two, rheostat controlled, 500K1V ^emersion heaters were 

used to heat the tank. Temperature control within 1/2°F. 
was provided by a Thermistemp unit controlling a third 
emersion heater rated at 750KW. The tank was kept at a 
uniform temperature by a variable speed electric stirrer. 

The tank was located on top of a second tank contain- 
ing the radioactive PU-BE source. This steel tank was 
2.95 feet in diameter, 1.25 high and contained 65 gals, of 
water. The sources were contained in three stainless steel 
tubes arranged in a triangular configuration.. The highest 
source was located 5/8 inches from the top of the tank. 

The sources were located directly under the foil column 
by three plexiglass extension rods. 

The three source strengths were 1, 1, and 3 curies. 
Since the sources were anisotropic, the extension rods were 
marked so that the sources could be placed to insure the 
same exposure during every run. Figure 2 shows the ex- 
perimental set-up. 
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Saturated boric acid solutions were added for poisoned 
runs. After thorough mixing, samples were drawn and an- 
alyzed cliecmically as described in Appendix A. 

Fifteen foils were used for each run. The relative 
flux was measured with bare indium foils 1.0 inches in 
diameter and .01 inches thick. Indium resonance activity 
was measured with similar foils placed in cadmium covers 
.020 inches thick. Irradiation time for each run was three 
hours. The foils were removed, one at a time, and allowed 
to decay for three minutes before counting. This waiting 
time allowed the 13 second half-life In^^^ to decay to a 
level that would not interfere with the 54 minute activit)'. 
Foils were counted for three minutes using the equipment 
specified in Table I. This counting technique resulted in 
standard deviations ranging from .1 to 10 percent. Back- 
ground counts were taken before and after each run and 
averaged. The average background was subtracted from the 
foil counts before correction factors were applied. 

Detector GM tube 



operating voltage 
Voltage Supply 
Preamplifier 



Hewlett-Packard Model 1045A 



Hewlett-Packard Model 6516A 



1300 volts 



Time-Scalar 



Hewlett-Packard Model 5201L 



TABLE I. 



COUNTING EQUIPMENT 
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FIGURE 2 EXPERIMENTAL APPARATUS 
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FIGURE 3. COUNTING EQUIPMENT 




IV. ANALYSIS OF DATA 



A. ACTIVATION DATA 

Each run was conducted with fifteen foils. The foils 
were irradiated for three hours, removed one at a time, 
then counted- starting with the foil closest to the source. 
In this way, the foils with the lowest counts were irradia 
ted as long as possible to improve the counting statistics 
Foil counts were corrected for differences in foil weights 
absorption of indium resonance neutrons by the cadmium and 
counter resolving time. 

Foil weights were normalized to 1.0 gram. Correction 
factors ranged from 1.071 to 1.082. 

The indium resonance counting rate correction was 
1.11 as determined from reference [10]. 

Counter resolving times were determined by the two- 
source method using 




1-N^t 1-N2t 



where A is the saturated activity, N is the mean counting 
rate, and t is the resolving time. A correction factor 
of 1/(1-N^) was applied to all counting times. 

The saturated activity was calculated using the cor- 
rected counts in the equation 

A = AC 

s ' 3-X(ta- t.) 
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where , 

ti - irradiation time 
tz - time counting started 
ts - time counting stopped 
tz~ti - \;aiting time 
t 3 ~tz - counting time 

Having calculated the relative saturated activity, 
the quantities Asr^ , As and A^^r"* were computed 

and plotted versus the distance from the source. The com 
puter program used for these calculations is described in 
Appendix D, 



B. RELAXATION LENGTH 

Equation (6) gives the space dependence of the slow- 
ing- down density as 

q (r,E) = — ^ 

Since the activity is proportional to the slowing 
down density as shown in Equation (21) , we- have 

C 

(32) 

Rewriting and taking the logarithum of both sides yields 



In A r^ 
s 

The quantity In A^r^ \>?as 
relaxation length is the 



= CONST - r/b 



(33) 



plotted against distance. The 
negative reciprocal of the slope 
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The same method was used to determine the relaxation 



length associated with the indium resonance activity. 



C. AGE AND MIGRATION AREA 

Equations (24) and (26) were used to determine Age 
and Migration Area. Computer Program No. 1 (Appendix D) 
plotted graphs of Asr^, A^^r^, Asr**, and A^^r"* versus r. 
A smooth curve was drawn through the data points up to a 
radius of 12 cm. At 12 cm the plot was considered expo- 
nential and extrapolated to infinity by means of a curve 
formed from the equation 



Asr^ = C e 



- r/b 



where b is the relaxation length, determined from curve 

fitting of measured activities and the logarithmic plot. 

The areas under the Asr^ and A ,r^ curves were determined 

cd 

by Simpsons rule. The areas for Asr** and A^^r"* were de- 
termined by multiplying the Ar^ data points by r^ and 
integrating the result, rather than integrating the Asr"* and 
A^^r"* curves directly. After the areas had been computed 
the second moments were determined by Equations (23) and 
(26) ; and finally Age and Migration Area were found from 
Equations (24) and (26) . The computer program used for 
these computations is Appendix E. Plots of Asr^ and A^^r^ 
are included in Appendix B. 



D. DIFFUSION LENGTH 

Diffusion length is determined from Equation (11) . 

The Age that has been determined is to the indium resonance. 
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1.8 cm^ was added to the Age to account for the Age in- 
crease from indium resonance to thermal energies [6]. The 
nature of the calculation- the difference of two relatively 
large numbers - leads to considerable variation in the re- 
sults even under the best experimental conditions. Small 
errors (±.5 cm^) in the computed values of the Migration 
Area and Age lead to 20 percent change in the value of 
diffusion length. If errors are estimated to be 2.5 cm^ , 
a very reasonable estimate for this experiment, the value 
of diffusion length at higher poison concentrations will 
vary 200 percent. As the poison concentration is increased, 
diffusion length decreases and using the same estimates of 
error the percent error becomes very large. From this it 
is concluded that the results of this experiment serve on]y 
to shoAv the trend of the temperature dependence of diffusion 
length and the values determined should not bo considered 
absolute . 
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TABLE II 



TEMPERATURE 95°F 



Boric Acid g/1 
Concentration 


Age 




Migration 

Area 


Diffusion 

Length 


8 ± .5 


58.5 ± 2.1 




62.5 ± 1.9 


1.5 + .6 

- 1.3 


16 ± .6 


58.6- ± 1.8 




61.6 ± 2.0 


1.1 + .9 

- 1.1 


30 ± 1.1 


58.4 ± 1.9 




60.9 ± 2.2 


.9 + 1.9 
- .9 




TABLE 


III 






TEMPERATURE 131° 


F 








Boric Acid g/1 
Concentration 


Age 




Migration 

Area 


Diffusion 

Length 


8 ± .5 


59.5 ± 1.8 




63.6 ± 1.7 


1.5+ .7 

- 1.5 


16 ± .6 


59.6 ± 1.9 




60.9 ±1.8 


1.2 + .6 
- 1.2 


30 ±1.1 


59.6 ± 2.0 




62.2 ± 2.2 . 


.9 + 1.1 
- .9 




TABLE 


IV 






TEMPERATURE 167° 


F 








Boric Acid g/1 
Concentration 


Age 




Migration 

Area 


Diffusion 

Length 


8 ± .5 


60.2 ± 1.7 




64.4 ± 2.1 


1.6 + .8 
- 1.6 


16 ± .6 


60.4 ± 1.7 




63.1 ± 2.0 


1.2+ .6 
- 1.2 



30 ± 1.1 62.9 ± 1.9 
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60.3 ± 2.2 



.9 + 1.0 
- .9 



and Tjj, versus POISON CONCENTRATION AT 95°F 
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POISON CONCENTRATION (g/I) 



and Ttm versus POISON CONCENTRATION AT 131°F 
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POISON CONCENTRATION (g/1) 
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POISON CONCENTRATION (g/1) 
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SQUARE DATA REF. [2] AT 22. 3° C 
CURVES FOR 22.3, 35, 55, 65, 70°C. 
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FIGURE 7 
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V. CONCLUSIONS 



Both Age and Migration Area were found to be relatively 
weak temperature dependent parameters in the temperature 
range under investigation. The trend of the temperature 
dependence was clearly established and in agreement with 
the theory. Errors associated with the calculations of 
these parameters were mainly due to the statistical errors 
from the weak neutron sources and the location of the neu- 
tron origin. The products A^r^, A^r“* , A^^r^, and A^^r“* 
were greatly influenced by the position of the origin. The 
triangular configuration of the sources made the selection 
of one origin for all neutrons difficult. The first origin 
selected was the geometric center of the three sources. 

The validity of this reference point vvras tested by com- 
paring the relaxation lengths obtained for pure water with 
those reported in references [6, 9]. The results were not 
in agreement. The best agreement was obtained by assuming 
an origin on the top surface of the uppermost (also the 
strongest) source. However, it should be pointed out that 
this forces the assumption that all neutrons are born at 
this point when actually the neutrons from the other 
source are born at a significantly different distance. 

The diffusion length in this experiment was determined 
by computing the square root of the difference between Age 
and Migration Area. Because of the nature of the calcula- 
tion the relatively small uncertainties associated with 
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these parameters can and did lead to very large uncertain 
ties in the diffusion length. The investigator is hesi- 
tant to place too much weight on the final values of dif- 
fusion length. Rather than draw any conclusions as to 
exact values of diffusion length the results are reported 
only to show the trend of the temjDerature dependence of 
diffusion length. 

Using a polynomial expansion suggested by Ref [2] a 
successful effort was made in representing both the temp- 
erature and the absorption dependence of diffusion length 
The form of the temperature dependence was derived from 
the MaxiArellian distribution and applied to the macroscopi 
absorption cross section. The plots of the approximation 
were in excellent agreement with the data reported in 
reference []]. It is hoped that this approach will serve 
as a basis for further work in this area and that the ap- 
proach will prove successful at higher temperatures and 
poison concentration. 
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APPENDIX A: DETERMINATION OF BORIC ACID CONCENTRATION 

A. 1 gram of Boric Acid (H^BO^) was dissolved in 300 ml 
of distilled H 2 O. Twenty-five ml of glycerol and 10 drops 
of phenophthlen indicator were then added and the result- 
ing solution was titrated using a .IN NaOH solution. A 
distinct reddish-pink color marked the titration end-point. 
This titration established the molarity of the standard. 

B. Sample Analysis 

A 25 ml sample was drawn for analysis. To this sam- 
ple were added 5ml of glycerol and 4 drops of phenophthalein. 
This sample was then titrated using a .IN NaOH solution. 

C. Sample Calculation 

143 ml of standard (.IN) NaOH solution was used to 
reach end-point of the standard solution (1 gm dissolved 
in 300 ml H 2 O) . 27.5 ml of NaOH produced end-point in 

sample. 



27.5 ml 

25 ml (SAMPLE) 



y ml 

300 ml (STANDARD) 



y = 33 ml 



From standard 



sol ' gm/liter 



143 ml 
3.33 g/1 



33ml 



X 



X = 7.7 gm/liter 
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APPENDIX B 



LOTS OF A and A FOR POISON RUNS 
s cd 
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APPENDIX C 



ACTIVATION DATA FOR POISON RUNS 
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TABLE V 



H,BO^ CONCENTRATION 8.0 g/1; INDIUM RESONANCE ACTIVATION 
TEMPERATURE 9 5°F 



r (cm) 


A^j^CCPM) 


MEAN COUNT 


STANDARD 

DEVIATION 


3.81 


4338.9062 


3510.6030 


59 . 25 


5.08 


2691.9985 


2172.0639 


46.60 


6.35 


2028.5357 


1631.6850 


40.39 


7.62 


1447.7858 


1160.9482 


34.07 


8.89 


1094.8874 


8713.8183 


29.51 


10.16 


769.0697 


612.2907 


24.74 


11.43 


580.6441 


460.4277 


21.45 


12.70 


423.2587 


-334.1813 


18.28 


13.97 


315.7761 


248.0522 


15.74 


15.24 


235.0387 


183.7874 


13.55 


16.51 


171.5650 


133.4304 


11.55 


17.78 


127.5265 


98.6032 


9.92 


19.05 


89 .6837 


68.9076 


8.30 


20.32 


67.0308 


51.1201 


7.14 


21.59 


75.2695 


40.5087 


6.36 



TABLE VI 

H BO, CONCENTRATION 8.0 g/1; TOTAL ACTIVATION 
TEMPERATURE 9 5°F 

STANDARD 



r (cm) 


Ag(CPM) 


MEAN COUNT 


DEVIATION 


3.81 


11509.8781 


13952.5742 


118.12 


5.08 


9497.4623 


8464.7929 


92.00 


6.35 


7423.1657 


6507.8710 


80.67 


7.62 


5566.3006 


4894.4101 


69.96 


8.89 


4071.2426 


3612.4257 


60.10 


10.16 


2972.1105 


2627.4343 


51.25 


11.43 


2250.4948 


1981.6005 


44.51 


12.70 


1729.0450 


1505.0419 


38.79 


13.97 


1235.7078 


1077.4638 


32.82 


15.24 


989.6018 


858.1154 


29.29 


16.51 


732.5951 


622.7670 


24.95 


17.78 


522.9899 


473.8107 


21.76 


19.05 


419 .0863 


356.7763 


18.88 


20.32 


322.2161 


272.4084 


16.50 


21 . 59 


248.4556 


208.4776 


14.43 
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TABLE VII 



H„BO„ CONCENTRATION 8.0 g/1; INDIUM RESONANCE ACTIVATION 
TEMPERATURE 131°F 



r ( cm) 


a^dCcpm) 


MEAN COUNT 


STANDARD 

DEVIATION 


3.81 


3740.3436 


3040.1044 


55.13 


5.08 


2271.1732 


1841.5061 


42.91 


6.35 


1736.7794 


1404.8835 


37.48 


7.62 


1290.1232 


1041.6748 


32.25 


8.89 


889.0577 


715.0063 


26.73 


10.16 


675.6499 


541.6247 


23.27 


11.43 


540.5715 


431.8361 


20.78 


12.70 


365.6290 


290.9912 


17.05 


13.97 


--269.6938 


^13.8565 


14.62 


15.24 


207.0582 


163.2117 


12.77 


16.51 


159 . 8396 


124.7460 


11.16 


17.78 


110.9198 


86.0412 


9.27 


19.05 


107.6441 


83.0443 


9.11 


20.32 


72.8171 


55.8133 


7.47 


21.59 


60.8034 


46.2799 


6 . 80 



TABLE VIII 



H„BO„ CONCENTRATION 8.0 g/1; TOTAL ACTIVATION 
TEMPERATURE 131°F 



r (cm) 


Ag(CPM) 


MEAN COUNT 


STANDARD 

DEVIATION 


3.81 


20109.0553 


18046.6484 


134.33 


5.08 


11566.5587 


10350.9726 


101.73 


6.35 


8869.6275 


7913.3593 


88.95 


7.62 


6536.7874 


5812.9609 


76.24 


8.89 


4962.5813 


4397.5390 


66.31 


10.16 


3720.1893 


3284.1132 


57.30 


11.43 


2812.3892 


2472.5913 


49 . 72 


12.70 


2091.6808 


1830.8740 


42.78 


13.97 


1610 . 7691 


1403.2451 


37.45 


15.24 


1271.6344 


1101.6093 


33.19 


16.51 


935.2484 


806.1274 


28.39 


17.78 


729.9680 


625.0952 


25.00 


19.05 


546.4106 


465.1699 


21.56 


20.32 


412.6101 


348.8295 


18.67 


21.59 


327.0331 


274.4113 


16.56 
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TABLE IX 



H^BO^ CONCENTRATION 8.0 g/1 ; INDIUM RESONANCE ACTIVATION 
TEMPERATURE 16 7° F 



r (cm) 


A^dCCPM) 


MEAN COUNT 


DEVIATION 


3.81 


4268.5135 


3436.8632 


58.62 


5.08 


2863.6486 


2299.7617 


47.95 


6.35 


2335.1988 


1870.2941 


43.24 


7.62 


1668.7725 


1332.7524 


36.50 


8.89 


1369.9945 


1090.5432 


33.02 


10.16 


973.7141 


772.9123 


27.80 


11.43 


697.1962 


551.6057 


23.48 


12.70 


475.3571 


374.7990 


19.35 


13.97 


388.4696 ' 


303.0793. 


17.40 


15.24 


392.8802 


255 .2706 


15.97 


16.51 


264.9628 


205.1085 


14.32 


17.76 


316.3616 


162.4721 


12.74 


19.05 


139.1435 


106.7854 


10.33 


20.32 


107.4450 


82.6591 


9.09 


21.59 


89.9828 


68.8904 


8.30 



TABLE X 



H BO^ CONCENTRATION 8.0 g/1; TOTAL ACTIVATION 
TEMPERATURE 16 7° F 



r ( cm) 


Ag (CPM) 


MEAN COUNT 


DEVIATION 


3.81 


20692. 8863 


18570.5976 


136.27 


5.08 


12481.0852 


11170.0156 


105.68 


6.35 


9169.2326 


8180.6640 


90.44 


7.62 


6277.5951 


5582.4726 


74.71 


8.89 


5344.6326 


4736.0898 


68.81 


10.16 


3937.0296 


3475.5358 


58.95 


11.43 


2886.7369 


2537.9560 


50.37 


12.70 


2290.5178 


2004.9187 


44.77 


13.97 


1768.2372 


1540.4257 


39.24 


15.24 


1368.2429 


1185.8757 


34.43 


16.51 


1108.0042 


955.0327 


30.90 


17.78 


861.6697 


738.2927 


27.17 


19.05 


659.6030 


561.5327 


23.69 


20.32 


477.2240 


403.4556 


20.08 


21.59 


346.5687 


290.9023 


17.05 
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TABLE XI 



H_BO CONCENTRATION 16 g/1 ; INDIUM RESONANCE ACTIVATION 
TEMPERATURE 95° F 



r(cm) 


AcdCcpm) 


MEAN COUNT 


STANDARD 

DEVIATION 


3.81 


4116.3652 


3151.4687 


56.13 


5.08 


3034.3682 


2334.3879 


48.31 


6.35 


2185.4791 


1689.1000 


41.09 


7.62 


1419.0948 


1101.6096 


33.19 


8.89 


1164.1740 


907.5075 


30.12 


10.16 


755.3704 


591.1831 


24.31 


11.43 


558.3584 


438.6540 


20.94 


12.70 


455.3512 


359.0251 


18.94 


-13.9 7 


309.0079 


244.4800 


15.63 


15.24 


227.6572 


180.7092 


13.44 


16.51 


220.7544 


175.7788 


13.25 


17.78 


209.2343 


167.1031 


12.92 


19.05 


131.9818 


105.7061 


10.28 


20.32 


130.0777 


104.4638 


10.22 


21.59 


117.1275 


94.3071 


9.71 



TABLE XII 



H BO_ CONCENTRATION 16 g/1; TOTAL ACTIVATION 
TEMPERATURE 9 5°F 



r (cm) 


Ag(CPM) 


MEAN COUNT 


STANDARD 

DEVIATION 


3.81 


18626.1055 


9497.1953 


97.45 


5.08 


10895.5391 


5582.4882 


74.71 


6.35 


6685.8054 


3441.4116 


58.66 


7.62 


5595.1975 


2890.9396 


53.76 


8.89 


3702.6471 


1922. 2897 


43.84 


10.16 


2939.2058 


1532.0256 


39.14 


11.43 


2357.4266 


1233.4499 


35.12 


12.70 


9306.7335 


804.2429 


28.35 


13.97 


1391.1851 


733.0488 


27.07 


15.24 


970.9109 


513.2778 


22.65 


16.51 


728.9893 


386.5915 


19.66 


17.78 


516.2110 


274.5700 


16.57 


19.05 


404.9721 


216.0157 


14.69 


20.32 


313.6666 


167.7666 


12.95 


21.59 


230.3158 


123.5047 


11.11 
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TABLE XIII 



H-BO^ CONCENTRATION 16 g/1; INDIUM RESONANCE 
TEMPERATURE 131°F 



r(cm) A(.p(CPM) COUNT 



3.81 


4133.0768 


3164.2629 


5.08 


2710.5077 


2085 .2368 


6.35 


2175.8190 


1681.6340 


7.62 


1594.8277 


1238.0266 


8.89 


1188.0543 


926 . 1230 


10.16 


876.8356 


686 . 2463 


11.43 


687.5132 


540.1201 


12.70 


501. 3977 


395 . 3308 


13.97 


382.6678 


302.7580 


15.24 


282.6519 


224.3627 


16.51 


224.8465 


179.0371 


17.78 


181.8621 


145.2425 


19.05 


127.9392 


102.4683 


20.32 


97.5775 


78.3633 


21.59 


73.5245 


59.1995 



TABLE XIV 



H-BO^ CONCENTRATION 16 g/1; TOTAL ACTIVATION 
TEMPERATURE 131°F 



r (cm) 


Ag (CPM) 


MEAN COUNT 


3.81 


15732.0113 


8021.5351 


5.08 


9167.6452 


4697.1757 


6.35 


6695.2128 


3446.2539 


7.62 


5151.0925 


2663.1140 


8.89 


3788.3850 


1966.8022 


10.16 


2832 . 8391 


1476.5834 


11.43 


2144.8511 


1122.2265 


12.70 


1595.8782 


838.0163 


13.97 


1177.5317 


620.4697 


15.24 


906.2015 


479.0686 


16.51 


668.0390 


354.2690 


17.78 


505.9825 


269.1296 


19.05 


354.5358 


189.1125 


20.32 


287.2258 


153.6246 


21.59 


213.6016 


114.5419 



ACTIVATION 



STANDARD 

DEVIATION 

56.25 

45.66 

41.00 

35.18 
30.43 

26 . 19 
23.24 
19.88 
17.39 
14.97 
13.38 
12.05 
10.12 

8.85 

7.69 



STANDARD 

DEVIATION 

89.56 
68.53 
5 8.70 
51.60 
44.34 
38.42 
33.49 
28.94 
24.90 
21.88 
18.82 
16.40 
13. 75 
12.39 
10.70 



58 



TABLE XV 



H,BO^ CONCENTRATION 16 g/1; 
TEMPERATURE 16 7° F 



r (cm) 


Acd(cpm) 


3.81 


4213.6345 


5.08 


2672.1145 


6.35 


2102.6545 


7.62 


1520 . 5024 


8.89 


1145.9158 


10.16 


869.8834 


11.43 


673.7225 


12.70 


492.3515 


13.97 


365 . 0609 


15.24 


278.6658 


16.51 


211.2086 


17.78 


158.9448 


19.05 


123.8967 


20 . 32 


90.8913 


21.59 


70 . 3152 



INDIUM RESONANCE ACTIVATION 

STANDARD 



MEAN COUNT 


DEVIATION 


3225.9375 


56.79 


2055 . 7004 


45.33 


1625.0871 


40.31 


1180 . 3295 


34.35 


893.2749 


29.88 


680 . 8051 


26.09 


529 . 2858 


23.00 


388.1984 


19.70 


288 . 8278 


16.90 


221.1986- 


14.87 


168.1763 


12.96 


126.9398 


11.26 


99.2307 


9.96 


72.9937 


8.54 


56.6154 


7.52 



TABLE XVI 



H,BO^ CONCENTRATION 16 g/1 
TEMPERATURE 16 7° F 



r ( cm) 


Ag(CPM) 


3.81 


15683.6585 


5.08 


9559.2448 


6.35 


7208.4950 


7.62 


5502.6384 


8.89 


3895.1716 


10.16 


2960.5430 


11.43 


2208.6044 


12.70 


1682.5217 


13.97 


1226 . 8919 


15.24 


936 . 2943 


16 . 51 


698.5095 


17.78 


500.3817 


19.05 


395.2241 


20.32 


299 . 4661 


21.59 


216.0595 



TOTAL ACTIVATION 

STANDARD 



MEAN COUNT 


DEVIATION 


7996.8828 


89.42 


4897.8164 


69.98 


3710.4577 


60.91 


2844.8632 


53.33 


2022.2421 


44.96 


1543.1474 


39.28 


1155.5834 


33.99 


883.5139 


29 . 72 


646.4787 


2-5.42 


494.9772 


22.24 


370.4277 


19.24 


266.1506 


16.31 


210 . 8160 


14.51 


160.1714 


12.65 


115 . 8599 


10.76 
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TABLE XVII 



H_BO_ CONCENTRATION 30 g/1; INDIUM RESONANCE ACTIVATION 
TEMPERATURE 9 5°F 



r (cm) 


A^p(CPM) 


MEAN COUNT 


STANDARD 

DEVIATION 


3.81 


3814 . 3995 


1752.1708 


41.85 


5.08 


2444.9645 


1128.5703 


33.59 


6.35 


1898.0324 


880 . 1640 


29.66 


7.62 


1466.2170 


682.9135 


26.13 


8.89 


1152.9376 


539.2492 


23.22 


10.16 


807.3318 


379.1101 


19.47 


11.43 


595.1468 


280.5332 


16.74 


12.70 


455.3512 


215.4151 


14.67 


13.97 


351.5190 


166. 8683 


12.91 


15.24 


238.8134 


113.7388 


10.66 


16.51 


181.2051 


86.5722 


9.30 


17.78 


125.2512 


60.0184 


7.74 


19.05 


88.8686 


42.7057 


6.53 


20.32 


44.0997 


21.2496 


4.60 


21.59 


17.3741 


8.3934 


2.89 



TABLE XVIII 



H BO- CONCENTIUTION 30 g/ 
TEMPERATURE 9 5°F 



r (cm) 


Ag(CPM) 


3. 


81 


9512. 


9576 


5. 


08 


5909 . 


7507 


6. 


3S 


4424. 


9625 


7. 


62 


3384. 


5104 


8. 


89 


2552. 


1833 


10. 


16 


1875. 


4493 


11. 


43 


1474. 


2346 


12. 


70 


1058. 


7 32 5 


13. 


97 


773. 


8019 


15 . 


24 


579 . 


9679 


16. 


51 


422. 


9913 


17. 


78 


322. 


3804 


19 . 


05 


242. 


7463 


20 . 


32 


202. 


5797 


21. 


59 


110. 


3097 



; TOTAL ACTIVATION 

STANDARD 



MEAN COUNT 


DEVIATION 


4850.5273 


69.64 


3027.9460 


55.02 


2277.6787 


47.72 


1749 . 7912 


41.83 


1325.0078 


36.40 


965 . 1225 


31.06 


767.2319 


27.69 


556.7866 


23.59 


408.3164 


20.20 


307.0202 


17.52 


225.1728 


15.00 


171.6843 


13.10 


129.6355 


11.38 


108.4727 


10.41 


59 . 2159 


7.69 



60 



1 



i 



t 



I 



I 













TABLE XIX 



H BO^ CONCENTRATION 30 g/1 ; INDIUM RESONANCE ACTIVATION 
TEMPERATURE 131°F 



r (cm) 


A^p(CPM) 


MEAN COUNT 


STANDARD 

DEVIATION 


3.81 


4032.5113 


3048.0051 


55 . 20 


5.08 


2261.6839 


2021.6328 


44.96 


6.35 


2034.2463 


1574.2250 


39.67 


7.62 


1441.6142 


1120.4519 


33.47 


8.89 


1119.1448 


871.3825 


29.51 


10.16 


845.8544 


660.1320 


25.69 


11.43 


639 . 5929 


501.1230 


22.38 


12.70 


483.9967 


380 .6340 


19.50 


13.97 


355.0934 


280.2563 


16.74 


15.24 


272.6566 


215.9255 


14.69 


16.51 


217.1437 


172.5205 


13.13 


17.78 


161.9829 


129.0929 


11.36 


19.05 


120.0964 


95.9932 


9 . 79 


20.32 


93.7458 


75.1415 


8.66 


21.59 


73.6595 


59.1995 


7.69 



TABLE XX 



H BO^ CONCENTRATION 30 g/1; TOTAL ACTIVATION 
TEMPERATURE 131°F 



r (cm) 


Ag(CPM) 


MEAN COUNT 


STANDARD 

DEVIATION 


3.81 


9348.0959 


4766.4648 


69.03 


5.08 


5167.0679 


2647.4218 


51.45 


6.35 


4059.7207 


2089.6765 


45.71 


7.62 


3025.5249 


1554 . 2170 


39.42 


8.89 


2255.1747 


1171.4912 


34.22 


10.16 


1693.3000 


883.1005 


29.71 


11.43 


1280.2822 


670.5717 


25.89 


12.70 


969.5395 


509.6274 


22.51 


13.97 


733.6620 


386.9528 


19.67 


15.24 


572.7375 


303.0563 


17.40 


16.51 


431.5252 


229.0411 


15.13 


17.78 


323.7266 


172.3308 


13.12 


19.05 


242.8408 


129.6355 


11.38 


20.32 


187.9917 


100.6241 


10.03 


21.59 


138.5682 


74.3591 


8.62 
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TABLE XXI 



H BO^ CONCENTRATION 30 g/1; 
TEMPERATURE 16 7°F 



r (cm) 


A^p(CPM) 


3.81 


3956.9440 


5.08 


2678.0214 


6. 35 


2029.5357 


7.62 


1547.6544 


8.89 


1174.0066 


10.16 


889.3505 


11.43 


684.7551 


12.70 


519.7670 . 


13.97 


393.3046 


15.24 


306.2453 


16 . 51 


239 . 2349 


17.78 


175.6441 


19.05 


137.7631 


20 . 32 


107.9603 


21.59 


80.4186 



INDIUM RESONANCE ACTIVATION 

STANDARD 



MEAN COUNT 


DEVIATION 


3029 . 4165 


55.04 


2060 . 2446 


45.38 


1568.5754 


39.60 


1201.4069 


34.66 


915.1723 


30.25 


696.0410 


26.38 


537.9533 


23.19 


409 . 8142 


20.24 


310.2597' 


17.61 


242 . 2943 


15.56 


189 . 8993 


13.78 


139 . 8590 


11.82 


110 .0232 


10.48 


85.8811 


9.26 


64.5832 


8.03 



TABLE XXII 



H-BO^ CONCENTRATION 30 g/1; 
TEMPERATURE 16 7° F 



r (cm) 


Ag(CPM) 


3.81 ■ 


8459.7042 


5.08 


5483.4342 


6.35 


4247.0050 


7.62 


3191.6215 


8.89 


2346.1385 


10.16 


1762.6128 


11.43 


1344.7866 


12.70 


997.1895 


13.97 


742.7238 


15.24 


588.4379 


16.51 


440.6199 


17.78 


91.4371 


19.05 


242.5287 


20.32 


200 . 2300 


21.59 


141.4392 



TOTAL ACTIVATION 

STANDARD 



MEAN COUNT 


DEVIATION 


4365.5546' 


66.07 


2841.7880 


53.30 


2209.9780 


47.01 


1656 .6054 


40.70 


1230 . 7382 


35.08 


927.8654 


30 . 46 


710.2744 


26.65 


528.3569 


22.98 


394.7204 


19.86 


313.6276 


17.70 


235.4891 


15.34 


48.9967 


6.99 


130 . 2845 


11.41 


107 . 8186 


10.38 


76 .3347 


8.73 
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APPENDIX D 



COMPUTER PROGRAM NO. 1 



THE PURPOSE OF THIS PROGRAM IS TO (L) CORRECT THE 
COUNTING DATA FOR NORMALIZED FOIL WEIGHT, ABSORPTION 
OF INDIUM RESONANCE NEUTRONS BY THE CADMIUM COVER AND 
COUNTER RESOLVING TIME, (2) COMPUTE THE MEAN COUNTING 
RATE AND, FROM THIS, DETERMINE THE STANDARD DEVIATION, 
(3) COMPUTE THE PRODUCTS A S*R**2 , A S=5=R**4 , LOG AS=!‘R')=»2, 
AS(CD)*R'i=-2, AS(CD)'-!=R«^4, AND (4) PLOT GRAPHS OF THE 
ABOVE PRODUCTS VERSUS R, THE DISTANCE FROM THE SOURCE. 



PROGRAM NOTATION 



D1 DECAY CONSTANT FOR INDIUM-116 

T COUNTER RESOLVING TIME 

POS DISTANCE FROM SOURCE 

X FOIL IRRADIATION TIME 

Y WAITING TIME 

Z WAITING TIME PLUS COUNTING TIME 

W CORRECTION FACTOR FOR FOIL WEIGHT 

Cl FOIL COUNTS 

SAT SATURATED ACTIVITY 

SIGMA STANDARD DEVIATION 

AMEAN MEAN COUNTING RATE 

F CADMIUM CORRECTION FACTOR 



IMPLICIT REAL*8( A-H,0-Z) 

WRITE (6,1) 

1 FORMAT(T5, 'DISTANCE' ,T25,* SAT ACT I V I TY ' , T5 6 , ' ME AN ' , 
1T87, 'SIGMA'//) 

REAL AC, AD,AA, BB,DISTA,DISTC, B,SUM,C,D, E,U , V, D1,C1, B1 
DIMENSION AC( 25) , AD( 25) , AA( 25) , BB( 25) , DISTA( 25) ,DISTC( 
25), B(3),SUM 

1(4),C(3,3),D(9),E(3) ,U(402 ) ,V(401) 

REAL LABEL1/4H /,LABEL2/4H /,LABEL3/4H /,LABE 
L4/4H / 

REAL48 ITITLE(12) 



RUN NUMBER THIRTEEN 



COMPUTATION OF SATURATED ACTIVITY OF INDIUM FOILS AND 
PRODUCTS AS*R=«=*2 AND AS4R**4 

Dl=1.28D-2 
T=6.0D-6 
D050 N=l, 15 

READ(5,2) POS,Cl,X,Y,Z 
2 FORMAT ( 5F9.0) 

READ (5, 6) W 
6 FORMAT ( 1F9.0) 

DIST=(POS-1.750)*0.254D01 
Cl = W*3.0*Cl/( 1.0D0-C1=!<T) 

ASAT=-Dl*Cl/( ( 1.0D0-DEXP(-D1*X) ) '4( D EXP ( -D1 *Y ) - DEXP ( -D 1« 
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Z) ) ) 

SAT=ASAT 
AMEAN=Cl/3.0 
SIGMA = AHEAN*^'0. 5D0 
F=1 . IIDO 

A=FSAT*DI ST=^--'2.0D0 

8 format ( 7?T12, F5l2?6X,D25.16,6X,D25. 16) 

WRITE (6f8) DIST,A,B1 
DISTC (N)=DIST 
AC(N)=A 
AD(N)=B1 
'50 CONTINUE 

rnMDiiTATTnM HF <^ATURATED ACTIVITY OF CADMIUM COVERED 
^8 'i[s 5nD PRODUCTS AS(C0)»R»»2 AND AS(CD)*P«»4 

D060 N = l,15 ^ V/ ^ 

READ(5,5) POStCI»X,Y,Z 
FORMAT ( 5F9. 01 
READ (5, 7) W 

01 ST= ( P0S-i?750 ) «0. 254D01 

ASAT*Di2c'r/((i:0D0-DE'xP(-Dl*X)l>MDeXP(-01*-Y)-DEXP(-Dl* 

Z) ) ) 

SAT=ASAT 
AMEAN=Cl/5.0 
9 I GMA=AMEAN**0 . 5D0 

WRITE (6t3) f>IST,SATtAMEAN,SIGMA 

FORMAT (T6tF 6.2, 6Xt D2 6. 16,6X,D25.16,6X,U26.io; 
A=SAT’!^DIST^!‘=i<2.0D0 
B1=SAT*DIST*-4.0D0 
WRITE(6,9) niST,A,Bl 

FORMAT ( /,T12, F6 .2,6X, D25.16,6X, D25. 16) 

DI STA( N)=DIST 
AA(N)=A 
BB(N)=B1 



5 

7 



PLOTTING OF COMPUTED PRODUCTS VERSUS R 



60 CONTINUE . . . r , , m 

READ (5,30) ( ITITLEI I ) ,1=1, 12) 

SlL^DRAW®(15,DISTC,AC,0,3,LABEL2,ITITLE,0,0,0,0,0,0,6 

REA0’*(5^3i) ( IT ITLE( I ) » 1 = 1 , 12 ) 

CAlT'^DRAW®(15,DISTC,AD,0,3,LABEL2,ITITLE,0,0,0,0,0,0,6 

REA0’('5^32) ( ITITLEd ) ,1=1,12) 

CALL^DRAW^ll5,DISTA,AA,0,4,LABEL3,ITITLE,0,0,0,0,0,0,6 

REAo’(5^33) ( ITITLE ( I ) , 1 = 1 , 12 ) 

CALrDRAW®(15,DISTA,BB,0,5,LABEL4,ITITLE,0,0,0,0,0,0,6 

.9,1, LAST) 

COMPUTATION/PLOTTING ^F LOG ^S*R^=*2 ^fjR a kp^cff ^ np * 
VERSUS R AND COMPUTATION OF A LEAST SQUARE . i 1 Ul- 

NOTE^^FOR^THE PURPOSE OF LEAST SQUARE FIT ONLY DATA 
POINTS FOUR THROUGH FIFTEEN ARE USED SINCE THIS IS 
THE ONLY PART OF THE CURVE USED TO ESTIMATE THE 
RELAXATION LENGTH 






D070 N=l,15 
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AC(iN)=ALQG(AC(N) ) 

AA(N)=ALOG( AA(NJ ) 

70 CONTINUE 

READ (5,34) ( IT ITLE (I ) , I = 1 , 12 ) 

34 F0RMAT(6A8) 

B( 1)=0.0 
3(2)=0.0 
3(3)=0.0 
D025 1=4,15 

B( 1) = B(1 )^•AC( I ) 

B(2)=B(2H AC(I)*0ISTC( I) 

B( 3) = 8( 3)+AC( I ) *DISTC( I )**2 

25 CONTINUE 
D026 L=l,4 
SUM (U =0.0 

26 CONTINUE 

DO 35 L=lr4 
DO 3 5 1=4, 15 

SUM(L) =SUM(L)+DISTC( I )**(L) 

35 CONTINUE 
C(1 V 1)=12 
D045 M=l,3 
D045 N=l,3 

IF(M+N.EQ.2) GO TO 45 
C(M,N)=SUM(M+N-2) 

45 CONTINUE 

WRITE(6,55> B( I ) ,B( 2 ) , B(3) 

55 F0RMAT(2X, 'VALUES OF B ' , lOX, 3F1 5,4} 

WRITE(6,65)(SUM(U ,L=l,4) 

65 FORMAT (2X SUM OF X I ' , lOX , 4F 1 5 . 4 ) 

WR I T E ( 6 , 7 5 ) ( ( C ( M , N ) , M= 1 , 3 ) , N= 1 , 3 ) 

75 F0RMAT(2X, 'VALUES OF C ' , 5X , 9F 12 . 4) 

CALL SIMQ(C,B,3,KS) 

WRITE(6,85) B( 1) ,B( 2) ,B(3) 

85 FORMAT (2X, 'SOLUTION' ,6X,3F10.4) 

U( 1 ) =10.0 
DO 95 1=1,401 

V( I ) =B( 1 )+U( I )=-'( B( 2)+U(I }*B( 3) ) 

U( I + 1)=U( I 1+0.05 
95 CONTINUE 

CALL DRAW (1 5 , D I STC , AC , 1 , 3 , LAB EL2 , I T I T LE , 0 , 0 , 0 , 0 , 0 , 0 , 9 
,15,1, LAST) 

CALL DRAW(401 ,U,V,2,0,LABEL2, ITITLE ,0,0, 0,0, 0,0, 6, 9,1, 
LAST) 

B( 1)=0.0 
B( 2 )=0.0 
B( 3 )=0.0 
D023 1=4,15 
B( 1 )=B( 1 )+AA( I ) 

B(2)=B(2)+AA( I )*DISTA( I ) 

B( 3 )=B( 3)+AA( I )*0I STA( I )’!==i'2 
23 CONTINUE 
D027 L=l,4 
SUM( L ) =0.0 

27 CONTINUE 
D03o L=l,4 
D036 1=4,15 

SUM(L) =SUM(L)+OISTA( I ) ( L ) 

36 CONTINUE 
C( 1 , 1) = 12 
D046 M=l,3 
D046 N=l,3 

IF(M+N.EQ.2) GO TO 46 
C(M,N)=SUM(M+N-2) 

46 CONTINUE 

WRITE(6,56) S( 1 ) ,B( 2 ) , B( 3) 

56 F0RMAT(2X, 'VZLUES OF B ' , lOX , 3F 1 5 . 4 ) 

WRITE(6,66MSUmL) ,L = 1,4) 

66 FORMAT ( 2X, ' SUM OF X I • , lOX , 4F 1 5 . 4 ) 

WRITE(6,76M ( C(M,N) ,M=1,3) ,N=1,3) 

76 F0RMAT(2X, 'VALUES OF C ' , 5X , 9F 1 2 . 4 ) 

CALL SIMQ(C,B,3,KS ) 
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WRITE(6t86) B( L) ,B(2) ,B(3) 

86 FORMAKZX, ‘SOLUTION' ,6X,3F10.4) 

U( 1 J=10.0 
D096 1=1,401 

V( I )=B( 1)+U( I )=!M B( 2)+U( I )*B(3) ) 

U( I + U =U( I )+0.05 
96 CONTINUE 

CALL DRAW(401 ,U, V,2,0 , LABE L3 , I T I TLE v 0 , 0 , 0 , 0 , 0 , 0 , 6 , 9 , 1 , 
LAST) 

CALL DRAW ( 1 5 , D I STA , AA , 3 ,4 , LABEL3 , 1 T I TLE , 0 , 0 , 0 , 0 , 0 , 0 , 6 
,9, 1,LAST) 

STOP 

END 
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COMPUTER PROGRAM NO. 2 



THE PURPOSE OF THIS PROGRAM IS TO COMPUTE DIFFUSION 
LENGtL THE AREAS UNDER THE CURVED PORTIONS OF THE 
CURVES AS=^R«*2,AS=i--R*-4.AS(CD)-R=i^=5=2 AND AS(CD)-R««4 ARE 
COMPUTED BY SIMPSON'S RULE. THE AREAS UNDER THE 
LINEAR PORTIONS OF THE CURVES ARE COMPUTED BY THE 
EQUATION Y-C2=!=EXP (-R/B ) . THE TOTAL AREAS ARE USED TO 
SOLVE FOR MIGRATION AREA AND AGE USING EQUATIONS (26) 
AND (24). DIFFUSION LENTH IS SOLVED BY 

DL= ( M I G- ( ACTE + l . 8 ) ) 5 WHERE THE 1.8 CORRECTION FACTOR 

CORRECTS FOR ADDITIONAL AGE REQUIRED TO REACH THERMAL 
ENERGIES. 



PROGRAM NOTATION 



B RELAXATION LENGTH ^ 

BCD RELAXATION LENGTH - CADMIUM COVERED FOILS 

R DISTANCE FROM SOURCE 

S2 AREA UNDER AS*R-«2 CURVE 

SC02 AREA UNDER AS ( C D ) -R'-^-2 CURVE 

S4 AREA UNDER AS*R-*4 CURVE 

SCD4 AREA UNDER AS ( CD ) -R*=«=4 CURVE 

MIG MIGRATION AREA 

AGE AGE 

DL DIFFUSION LENGTH ^ 

REAL*8 H, SUM2, SUM4,SCD2,S4, SCD4, A2 » ACD2 t AS 2 , AS4» ASCD2 
XS2 Y YCD C2tR 

DIMENSION Y{ 100) ,YCD( 100) ,R( 100) 



RUN NUMBER SIX - CONC 8 GM/L - TEMP 131F 



READ(5»10) ( Y( I ) »I = 1,10I 
10 FORHAT( IFIO.O) 

B=11.4 

C2=3000000.0 

READ(5,50) (YCD( I ) , 1=1 ,10) 
50 FORMAT! IFIO.O) 

BCD=9.0 

CD2=710000.0 



COMPUTING SATURATED ACTIVITY TIMES DISTANCE SQUARED AR 
AS* R^’^ 2 



D0I60 1=11,80 
R( I ) =1 .0*( I ) 

Y( I )=C2*DEXP(-R( I )/B) 

H=1.0 

N=79 

SUM2=0.0D0 

SUM4=0.0D0 

M=N^1 

160 CONTINUE 

D020 1=3, N, 2 
SUM2=SUM2+Y( I ) 



ooo ooo non ooooo ooo 



20 CONTINUE 

D030 1=2, H»2 
SUM4=SUM4+Y( I ) 

30 CONTINUE 

S2 = H/.3D1=5=(Y( I ) + .4Dl’!=SUM4+.2Dl=t<SUM2 + Y(N+l) ) 

WRITE(6,40) S2 

40 F0RMAT(T19, • SI MPSONS A PPROX I MAT I ON • , 025 . 16/ // ) 

AS(CD)*R=«'-=S'2 

D0230 1=11,80 
R( I )=1.0=J=( I) 

YC0{ I )=CD2*DEXP(-R( I )/BCD) 

H=1.0 

N=79 

SUM2=0.0D0 

SUM4=0.0D0 

M=N-1 

230 CONTINUE 

0060 1 = 3, N, 2 
SUM2=SUM2+YC0{ I ) 

60 CONTINUE 

0070 1=2, M, 2 
SUM4=SUM4+YC0( I ) 

70 CONTINUE 

SC02 = H/.301^= (YCD( I ) + . 40 1*SUH4+ . 201*SUM 2+YCO ( N + 1 ) ) 

WRITE (6,80) SCO?. 

80 F0RMAT(T19, • SIMPSONS APPROX I MAT I ON », 02 5 . 16/// I 

COMPUTING SATURATED ACTIVITY DISTANCE TO THE FOURTH AR 

AS*R«-4 

0090 1=1,10 
R( I ) = 1.0^M I) 

Y( n=Y{ I)=:R( I )*=«=2.0 
90 CONTINUE 
H=1.0 

EXPODENTIAL PORTION 

DO 200 1=10,80 
R(I) = 1.0=!'(n 

Y( n=C 2 *DEXP(-R( D/Bl’i'Rl I )=;=«2 .0 
200 CONTINUE 
N=79 

SUM2=0.0D0 

SUM4=0.0D0 

M=N-1 

00100 1=3, N, 2 
SUM2=SUM2+Y( I ) 

100 CONTINUE 

00110 1=2, M, 2 
SUM4=SUM4+Y{ I ) 

110 CONTINUE 

S4=H/e3Dl=<‘(Y( I ) + .4Dl=5^SUM4+.2Dl^SUM2+Y( N+U ) 
WRIT£(6,120) S4 

120 F0RMATIT19, ‘ SIMPSONS APPROX IMATION' , 025. 16/// I 

AS(CD)*R**4 

H=1.0 

N=15 

DO 130 1=1,10 
R{ I )=1 .0«( I ) 

YCD( I )=YCD( I )*R{ I )**2.0 
130 CONTINUE 

EXPODENTIAL PORTION 

DO 260 1=10,80 

R( n=i . 0 «( I ) 
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YCD( I )^C02*DEXP(-R( I )/BCD)*R( I )«=«'2.0 
260 CONTINUE 
N=79 

SUM2=0„0D0 

SUM4=0,0D0 

M=N-1 

D0140 1=3, N, 2 
SUM2=SUM2+YCD( I ) 

140 CONTINUE 

D0145 1=2, M, 2 
SUM4=SUM4+YCD( I ) 

145 CONTINUE 

SCD4=H/.3D1=!=(YC0( I ) + . 401=5=5 UM4+ . 2D1 * SUM2+YCD ( N + 1 ) ) 
WRITE(6,150I SCD4 

150 FORMAT (T19, • SI HPSONS A PPROX I HAT I ON • , 02 5 . 16/ / / J 



COMPUTING OIFFUSION LENGTH 

AS2=S2 

AS4=S4 

ASCD2=SC02 

ASC04=SC04 

MIG=AS4/( AS2=!=6.0) 

WRITE! 6,2901 MIG 

290 FORMAT (T19, * MIGRAT ION ARE A= ' , 02 5 . 1 6/ // ) 
AGE = ASC04/ (ASCD2=!-'6. OJ 
WRITE(6,300! AGE 

300 FORMAT (T 19, *AGE=* , 025.16/// » 

DL=(MIG-(AGE + 1.8) l*=«=.5 
WRITE(6,310) OL 

310 F0RMAT(T19, 'OIFFUSION LENGTH= * , 025 . 16/ // ) 
STOP 
ENO 
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APPENDIX F 



COMPUTER PROGRAM NO. 3 



THE PURPOSE OF THIS PROGRAM IS TO COMPUTE AND PLOT 
A POLYNOMIAL APPROXIMATION OF 1/L*=<'2 IK**2) AS A 
FUNCTION OF BORIC ACID CONCENTR ATI C3N AND TEMPERATURE. 
THE BORIC ACID CONCENTRATION IS EXPRESSED AS THE 
MACROSCOPIC ABSORPTION CROSS-SECTION OF BORON. THE 
TEMPERATURE DEPENDENCE FACTOR ( TN/ TO ) . 5 IS BASED ON 
A REFERENCE TEMPERATURE OF 22.3 DEGREES CENTIGRADE. 
REFERENCE DATA AT THIS TEMPERATURE IS FROM REF. (2). 



PROGAM NOTATION 



EA MACROSCOPIC ABSORPTIN CROSS-SECTION 

DL DIFFUSION LENGTH 

TN MODERATOR TEMPERATURE 

TO MODERATOR REFERENCE TEMPERATURE 



REAL DL,EA,B,C,U,Vf SUM,TN, TO,A 
REAL LABEL2/AH / 

DIMENSION DL(IO) ,EA(10) ,B(4) ,SUM(4) , U { 402 ) » V ( 40 1 ) , C ( 4 , 
4) , A(4) 

REALMS ITITLE (12) 

RE AD (5, 10) ({DL{I),EA(I)i,I=l,10) 

10 FORMAT ( 2F9.0) 

READ(5,34) ( IT I TL E ( I ) , I = 1 , 1 2 ) 

34 F0RMATI6A8) 

B(1)=0.0 
B(2)=0.0 
B(3)=0.0 
B(4)=0 .0 
D025 1=1,10 
DLI I ) = 1 .0/DL( I ) 

B(2) = B(2)+DL(n*EA( I ) 

B(3) = B{3)+DL( n^J'EAI I )^'^2 
B(4)=B(4)+DL{ I )=i=EA( I ) **3 

25 CONTINUE 
D026 L=l,6 
SUM(L)=0.0 

26 CONTINUE 
D030 L=l,6 
D030 I=i,10 

SUM! L) =SUH(U+EA( I )=J==^=(L) 

30 CONTINUE 
C(l, 1)=10 
DO 40 M=l,4 
DO 40 N=l,4 
IFIM+N.EQ.2) GO TO 40 
C(M, N)=SUM(M+N-2) 

40 CONTINUE 

WRITE{6,50) B( 1 ) , B( 2) , B(3) , B(4) 

50 F0RMATI2X, • VALUES OF B * , lOX , 4F 1 5 .4 ) 

WRITE(6,60) (SUM(L) ,L=1,6) 

60 FORMAT (2X, • SUM OF X I ' , lOX , 6F 1 5 . 4 } 

WRITE(6,70)( (C{ M,N) ,M=1,4) rN=l,4) 

70 F0RMATI2X, 'VALUES OF C',16F12.4) 
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CALL SIMQ(C,B,4,KS) 

WRITE(6,80) B( 1 ) ,B(2 ) ,B(3) , B(4) 

80 F0RMAT(2X, 'SOLUTION' ,6X»4F10.4) 

U( 1)=.01 
DO 90 1=1,401 

V( I )=B{ 1 H-U( I )^M B( 2)+U( I )=:=B( 3 ) + ( U( I )*='-2) *B( 4) ) 

U( I+1)=U( I )+.001 

90 CONTINUE 

CALL DRAW( 10 , E A, DL , 1 , 3 , LAB EL2 , I T IT L E , 0 , 0 , 0 , 0 , 0 , 0 , 6 , 9, 1 
,LAST) 

CALL DRAW(401,U, V,2,0,LABEL2,ITITLE,0,0,0,0,0,0,6,9, 1, 
LAST) 

TN=308.0 
T0=296.15 
A( 1 ) = B( 1 ) 

A(2)=B(2) 

A(3)=B(3) 

A(4)=B{4) 

A(4)=( (T0/lN)**.5)*A(4)*c8862 
A(3) = ( ( TO/TNI^i^’f^.B) -A( 3 8862 

A ( 2 ) = ( ( TO/ TN ) *=■:-' . 5 ) *A ( 2 ) * . 8 862 
A(l)=< (T0/TN)**.5)*A(l)'-!'.8862 
DO 91 1=1,401 

V( I )=A( i)+U( I )*( A( 2) +U( I )'^A(3)+( U{ I ) =:<*2 ) ^^'A ( 4 ) ) 

U( I+lf =U( I )+.001 

91 CONTINUE 

CALL DRAV>I(401 , U , V , 2 , 0 , LAB EL 2 , I T I TL L , 0 , 0 , 0 , 0 , 0 , 0 , 6 , 9 , 1, 
LAST) 

TN=328.0 

A(1)=B(1) 

A( 2) = B( 2) 

A(3)=B(3) 

A(4)=0(4) 

A(4) = { ( TO/TN)=!^*. 5)’}=A(4)*.8 862 
A<3) = ( ( T0/TN)4'^. 5) =?=A( 3 )4-. 8862 
A ( 2 ) = ( ( TO/TN ) . 5 ) =4=A ( 2 ) . a 862 

A(l) = ( (T0/TN)*=:^.5)*A( 1)*.8862 
DO 92 1=1,401 

V( I )=A( 1)+U( I )*( A( 2)+U( I )*A(3) + (U{ I ) **2 ) >.’=A ( 4 ) ) 

U( I+1)=U( 1 )+.001 

92 CONTINUE 

CALL DRAW(401 ,U,V,2,0, LABE L2, 1 TITLE, 0,0,0, 0,0, 0,6,9, 1, 
LAST ) 

TN=338.0 
A( 1)=B( 1) 

A(2)=B(2) 

A(3)=B(3) 

A(4)=B(4) 

AI1) = ( (T0/TN)**.5)*A( 1 )*.8862 
A<2) = { (T0/TN)=^*.5)=!^A(2)*.8862 
A(3) = ( (T0/TN)**.5)*A(3)=!=. 8 862 
A(4) = { (T0/TN)*={=.5)»A(4)*.8 862 
DO 93 1=1,401 

V(I)=A(1)+U(I)=^(A(2)+U(I)*A(3) + (U( I ) =«^*2 ) *A { 4 ) ) 

U( I + l ) =U( I ) + .001 

93 CONTINUE 

CALL DRAW(401 ,U,V, 2,0, LABE L2, 1 TITLE, 0,0, 0,0,0, 0,6,9, 1, 
LAST) 

TN=343.0 
A( 1)=B( 1) 

A(2)=B(2) 

A( 3)=B( 3) 

A(4)=B(4) 

A( 1 ) = ( (T0/TN)=^^*5)*A( 1 )*.8 862 
A(2)=( (TO/TN)**.5)*A(2)=5'-.8862 
A(3) = ( (T0/TN)*=!=.5)=!=A{ 3)*.8862 
A(4) = { (T0/TN)=i'*.5 )*A(4)*.8 862 
DO 94 1=1,401 

V( I )=A( 1 1-'-UI I A( 2) +U( I )*A(3 ) + ( U( I )4'*2) *A(4) ) 

U( I + 1)=U( I ) + .001 

94 CONTINUE 

CALL DRAW(401 ,U, V,3,0,LABEL2, ITITLE, 0,0,0,0,0,0,6,9, 1, 
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STOP 

END 
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